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We study Fermi edge singularities in photo-absorption spectra of generic mesoscopic systems such as quan- 
tum dots or nanoparticles. We predict deviations from macroscopic-metallic behavior and propose experimental 
setups for the observation of these effects. The theory is based on the model of a localized, or rank one, pertur- 
bation caused by the (core) hole left behind after the photo-excitation of an electron into the conduction band. 
The photo-absorption spectra result from the competition between two many-body responses, Anderson's or- 
thogonality catastrophe and the Mahan-Nozieres-DeDominicis contribution. Both mechanisms depend on the 
system size through the number of particles and, more importantly, fluctuations produced by the coherence char- 
acteristic of mesoscopic samples. The latter lead to a modification of the dipole matrix element and trigger one 
of our key results: a rounded K-edge typically found in metals will turn into a (slightly) peaked edge on average 
in the mesoscopic regime. We consider in detail the effect of the "bound state" produced by the core hole. 

PACS numbers: 73.21. -b,78.70.Dm,05.45.Mt,78.67.-n 



I. INTRODUCTION 

Fermi edge singularities (FES) are among the simplest 
many-body effects found in condensed matter physics, and 
have been studied extensively for bulk systems.^* 2 ^ The con- 
tinuing progress in the fabrication and experimental investiga- 
tion of mesoscopic systems makes probable that such singu- 
larities could be observable relatively soon in quantum dots or 
nanoparticles, for which interference effects, and thus meso- 
scopic fluctuations, have to be taken into account. The study 
of Fermi edge singularities in this mesoscopic regime is the 
topic of the present paper; it extends and deepens our previ- 
ous results 4£ 

For macroscopic systems, the physics we discuss is known 
as the "x-ray edge problem" and is well established and 
understood. 12,3 In an x-ray absorption process, a core elec- 
tron is excited into the conduction band. It leaves behind a 
positively charged hole that can be considered as a static im- 
purity. For electronic densities corresponding to an interaction 
parameter r s ~ 1, as realized typically in both semiconductors 
and metals, the screening length is of order the Fermi wave- 
length. This impurity can therefore be considered as local- 
ized (point like) once screening by the conduction electrons 
is taken into account. A similar situation arises after the ex- 
citation of a valence electron in semiconductor photolumines- 
cence studies in which recombination occurs at an impurity. 

The conduction electrons respond to such an abrupt, non- 
adiabatic perturbation by slightly adjusting their single parti- 
cle energies and wave-functions. Although the overlap of the 
single particle states before and after the perturbation is very 
close to one, the overlap A between the initial and final many- 
body ground states tends to zero in the thermodynamic limit 
as a power of the number M of particles. This effect is known 
as the Anderson orthogonality catastrophe (AOC). 6 As a re- 
sult of AOC, the photo-absorption cross section A(ui) will be 
power-law suppressed for energies to near the threshold en- 



ergy (Fermi edge) uj th . 

In the x-ray edge problem, AOC competes with a second, 
counteracting many-body response, also related to screening 
of the impurity potential by the conduction electrons. It is of- 
ten referred to as Mahan's exciton, Mahan's enhancement, or 
the Mahan-Nozieres-DeDominicis (MND) contributionji 2 - In 
contrast to AOC which acts universally, the MND response 
depends on fulfilling dipole selection rules. It therefore de- 
pends on the symmetry of the conduction and core states. We 
will distinguish two symmetries of the core electron wave- 
function throughout the paper: s-like symmetry (correspond- 
ing to the K-shell of the atom) and p-like symmetry (£2,3- 
shell). The respective thresholds are known as the K- and L- 
edge. 

The bulk x-ray edge problem was analyzed in detail us- 
ing different techniques. Early approaches by Mahani and 
Nozieres and co-workers 2 treated it based on field-theoretical 
methods and diagrammatic perturbation theory in analogy 
to the Kondo problem. 7 Schotte and Schotte 8 used a Fermi 
golden rule approach and employed bosonization techniques 
for the rotationally invariant, effectively one-dimensional, 
scattering potential. The x-ray edge problem was also ad- 
dressed by, among others, Friedel 9 and Hopfieldfi. In the 
1980's Tanabe and Ohtaka demonstrated in detail the suitabil- 
ity of a Fermi golden rule approach. 311 This method explic- 
itly uses the fact that the impurity in the x-ray edge problem 
is static rather than dynamic, as in the Kondo problem. The 
spirit of the Fermi golden rule approach is illustrated schemat- 
ically in Fig. [lj it is the approach used here. 

The effect of a local perturbation such as the one induced by 
the core hole is characterized by the partial-wave phase shifts 
Si (at the Fermi energy) for each orbital channel. The phase 
shifts have to obey the Friedel sum rule Z = J2i 2(2/ + l)5i/ir 
with the screening charge Z = — 1 in our case. The factor of 
two accounts for spin. 

The result of all the bulk techniques is that the photo- 
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absorption cross section A(ui) near threshold oj t h has the form 

ai \ i x-2^+y 2 (2z+i) ri-i 2 

The first term in the exponent involves only the optically 
active channel (labeled Iq) and is the MND contribution, 
whereas the second term sums over all channels and corre- 
sponds to AOC. Note the different functional dependence on 
the phase shifts, linear and quadratic, respectively. 

We will, throughout this paper, assume that (i) the local 
part of the conduction electron's wave-function - at the lattice 
level - is featureless (i.e. of s-type) and (ii) the perturbation 
created by the core hole is spherically symmetric. In the bulk, 
it follows from the first assumption that the optically active 
channel, for which core and conduction electrons are linked 
by the dipole operator, is Iq — 1 for K-shell core electrons and 
lo = for the L-shell. As a consequence, for the K-shell, the 
perturbation acts on a channel which is not optically active. 
Therefore, the absorption spectrum is only affected by AOC in 
the I — channel, yielding a suppression, or rounding, of the 
edge of the photo-absorption spectra ("rounded edge"). On 
the other hand, for L-shell core electrons, the optically active 
channel is the one affected by the perturbation, and the cor- 
responding materials typically show an enhancement of the 
photo-absorption at the threshold frequency ("peaked edge"). 
For a detailed analysis of the x-ray spectra of bulk Li, Na, Mg, 
and Al, we refer the reader to Ref. 12, where in addition ef- 
fects due to phonon excitations, the finite lifetime of the hole, 
and the deviation from spherical symmetry of both the local 
part of the conduction electron state and the perturbing poten- 
tial are considered, all of which we will neglect here in order 
to focus on the essential mesoscopic physics. 

In semiconductors, typically only s-like conduction elec- 
trons exist which, consequently, have to provide all the screen- 
ing of the core hole. From the Friedel sum rule we then find 
5q = — it/2. This implies that we are in the strong perturbation 
regime which is typically not realized in bulk metals or related 
nanoscale structures like metallic nanoparticles. We shall see 
below that this has significant physical consequences related 
to the formation of a bound state. 

In contrast, in metals, electrons of all channels Si typi- 
cally contribute to Friedel screening. Thus the phase shift 
in the optically active channel is not tt/2. However, the x- 
ray edge physics can be successfully captured based on 5q 
alone, i.e., by assuming a spherically symmetric core hole 
potentia h 2 ' 3 ' 12 ' 13 This remains true in the mesoscopic regime. 
For this situation, then, one considers do < 7r/2 even though 
there is only one phase shift in the model. Thus, to address 
both the semiconductor and metal situations, we include re- 
sults for the full range <5o < tt/2. 

The above description applies to clean bulk systems. The 
question we would like to address in this paper is whether 
the results found for the edge behavior in metals also hold, or 
have to be modified, when considering smaller, fully coherent, 
mesoscopic or even nanoscopic samples. Here, we study the 
universal class of chaotic ballistic systems, to which a random 
matrix model of the energy levels and wave-functions applies. 

Reducing the size of the system will affect the Fermi edge 
singularity in various ways. First of all, the number of par- 
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FIG. 1: (Color online) Schematic illustrating processes contribut- 
ing to the photo-absorption cross section in a Fermi golden rule 
approach. Conduction electrons (mean level spacing d) in a 
generic chaotic system initially occupy levels eo • • • £m-i (filled) 
and em • • • £iv-i (empty), that lower to Ao . . . Aat-i when the core 
electron c is excited. Optically active electrons contribute via the co- 
herent superposition of (a) direct and (b) replacement processes. In 
addition, one (or more) electron-hole pairs can be generated in shake- 
up processes, (c) for optically active channel and (d) for spectator 
channel, which are especially important away from threshold. Their 
presence reflects the suddenness of the perturbation. The vertical ar- 
row in (a) is the bare process which represents the photo-absorption 
cross section in the naive picture without many-body effects. It de- 
pends only on the dipole matrix element w C j between the core elec- 
tron c and the single particle state j, A b (uj) oc \w c j | 2 . The effect of 
AOC is accounted for in the direct process by the additional factor 
|A| 2 < 1, A d (uj) oc |w C7 | 2 |A| 2 , 



tides (electrons) M in the system will be finite - we are not 
in the thermodynamic limit any longer - and for instance the 
power law dependence of A(M) causes, in fact, a huge dif- 
ference in the efficiency of AOC in systems with M ~ 10 23 
from those with, say, M ~ 100 electrons. Secondly, as the 
system becomes fully coherent, a plane wave (Bloch wave) 
description of the conduction electrons does not apply any 
longer. Rather, we need the actual wave-function of the spe- 
cific mesoscopic systems to be described. One consequence 
is that, since the confining potential will in general destroy 
spherical symmetry, angular momentum is lost as a quantum 
number, and the dipole selections rules need to be modified. 
Furthermore, interference effects, and therefore mesoscopic 
fluctuations, need to be taken into account. All of these affect 
both the Anderson orthogonality catastrophe and the Mahan- 
Nozieres-DeDominicis contribution. 

Aspects of AOC in disordered mesoscopic systems have 
been addressed in Refs.[l4j and [lj. AOC and x-ray photoe- 
mission spectra (where the excited core electron leaves the 
metal and the edge behavior is determined by the AOC re- 
sponse alone) were studied in Ref.[l6| for impure simple met- 
als. AOC in ballistic chaotic systems was the subject of the 
first paper 5 in this series. Our main findings - in line with the 
results in Refs.fTil and[l5l - are (i) incompleteness of AOC due 
to the finite number of particles, and (ii) a broad distribution 
of AOC-overlaps as a result of mesoscopic fluctuations which 
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(iii) are dominated by the levels around the Fermi energy. 

For (mesoscopic) x-ray absorption or photoluminescence 
spectra both the AOC and MND effects are of importance. 
The MND response will, of course, also depend on the sys- 
tem size. However, the relative strength of the two com- 
peting many-body responses might change as the system 
shrinks. Optical FES in one-dimensional quantum wire sys- 
tems have been studied both in experiment 17 and theory. 18 
In contrast, FES in the photo-absorption spectra of two- or 
three-dimensional mesoscopic system have, to the best of our 
knowledge, not been addressed in literature. 

The present paper aims at filling this gap. It is organized 
as follows. In Sec. [Tj] we introduce our model for the (rank 
one) perturbation of chaotic conduction electrons described 
by random matrix theory, the dipole matrix element, and how 
spin is taken into account. In Sec. [HI] we consider in more 
detail the formation and role of the bound state appearing for 
strong perturbations. In Sec.[IV]we explain our (Fermi golden 
rule) approach to the photo-absorption cross section. The re- 
sults are presented in Sec.|VI]for the average photo-absorption 
cross sections at the K- and L-edges, and in Sec. IVIII for the 
mesoscopic fluctuations. We devote Sec. I Villi to the discus- 
sion of feasible experimental setups that would allow probing 
of our results, and close with a summary in Sec.lIXI 



II. MODEL 

A. Initial and Final Hamiltonian 

In our model of a quantum dot or nanoparticle, the elec- 
trons are confined in a coherent, irregularly shaped system. 
We describe the unperturbed system (the conduction electrons 
before a core electron is excited) by the Hamiltonian 



which we consider here (as realized typically in both semicon- 
ductors and metals), the rank one approximation is, however, 
appropriate. The strength of the interaction between the core 
hole and the conduction electrons is quantified by the param- 
eter v c < 0, which in turn is related to the phase shift in the 
band center, So, by^ 3 . 



do = arctan — — 
a 



(3) 



In addition to v c , the effectiveness of the perturbation depends 
also on the wave-functions' amplitude (fik(f c ) at the position 
of the perturbation, which, in a mesoscopic system, will vary 
from state to state. In terms of Uk = \/Tltpk{r c ) so that 
(|ufc| 2 ) = 1 (with ft denoting the volume in which the elec- 
trons are confined), the perturbation can be expressed as 



■ 2_^, u k u k>C' k C k > 
kk' 



(4) 



For reasons of comparison, we also define the bulk-like situ- 
ation: equidistant unperturbed energy levels a distance d apart 
and constant uj. = 1 throughout the sample. More details 
concerning this model are given in Sec. II of Ref. 5 (the first 
paper of the series). 

Introducing cjl a as creator of a particle in the perturbed 
orbital ip K (f), we obtain the diagonal form of the perturbed 
Hamiltonian 



H 



(5) 



■ Cfc. 



(2) 



(Note that we will use Greek letters to refer to the perturbed 
system.) In analogy to Uk, we refer to the perturbed ampli- 
tudes by u K = \/VLip K (f c ). Finally, we define the transforma- 
tion matrix a = {a K k), 



with discrete eigen-energies e& (k = 0, . . . , N — 1). The op- 
erator ct a creates a particle with spin a =f , J, in the orbital 
ipk (r) ■ The energy levels of the unperturbed system follow the 
statistics of random matrix theor y 19,20 (RMT) and are charac- 
terized by a mean level spacing d, cf. Fig. Q] We will dis- 
tinguish situations where time-reversal symmetry is present 
(circular orthogonal ensemble, COE) from those where it is 
broken by, e.g., the presence of a magnetic field (circular uni- 
tary ensemble, CUE). 19,20 As the number of electrons does not 
change in the processes that we consider, we drop the charging 
energy term normally present in isolated mesoscopic systems. 
Furthermore, we neglect any change in the residual electron- 
electron interactions. 21 

This initial situation is perturbed by a rank one or contact 
potential V c acting at the location r c of the core electroni 2 ^ 
Models for perturbations which are more general than rank 
one can also be considered. 22 It is necessary to consider them, 
for instance, for high density electron gas (r s <C 1) for which 
the screening length is significantly larger than the Fermi 
wavelength. For the density range corresponding to r s ~ 1 



N-l 

^ a K k^Pk 

k=0 



(6) 



for later use. For relations between Uk, u K , ak K , and the e^, A K 
seeEqs. (15)-(19) in Ref.H 

We denote the M -particle, Slater-determinant ground states 
of the unperturbed and the perturbed system by |$o) ar, d |^o)> 
respectively. Their overlap is the Anderson overlap A which 
we considered in detail in Ref .|H In the case of a rank one per- 
turbation, it can be expressed as a function of the unperturbed 
and perturbed energy levels alone, 3 



M-l JV-1 



iAi 2 =K* i<i> )i 2 = ii n 



i=0 j=M 



(Aj -£j)(ej -Xj) 
(A, - Xi)(ej - a) 



(7) 



Note that whenever possible, we use the index j for levels 
above Ep (7 for perturbed levels), and i (p) for levels below 
Ep. Furthermore the index k (k) is reserved for reference to 
all unperturbed (perturbed) levels. 
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B. Dipole Matrix Element 

The dipole operator is D = (eEo/cm)(e- p + h.c.) where 
E is the magnitude of the electric field and e is its polariza- 
tion. The dipole matrix element 



(8) 



is therefore proportional to the overlap of the perturbed wave- 
function i/jp with the derivative along e of the core electron 
wave-function \<p c ). For K-shell core electrons (spherically 
symmetric), ip c is, on the scale of the Fermi wavelength of the 
conduction electrons, essentially a ^-function. In that case, 
w cfl is proportional to the derivative of along e. For L- 
shell core electrons on the other hand, the derivative of \<p c ) is 
approximately a <5-function, and therefore w cf2 is in this case 
proportional to ip^ itself rather than to its derivative. As a 
consequence, one has 



K-edge 
L-edge . 



(9) 



In the bulk, this implies precisely the selection rules men- 
tioned in the introduction, namely that the optically active 
channel (with non-zero w cfl ) is I = 1 for K-shell core elec- 
trons and I = for L-shell core electrons. Since only the I = 
channel is affected by the rank-one perturbation in Eq. ©, 
there is no MND enhancement for the K-shell, and thus it has 
a rounded Fermi edge due to AOC. 

In the mesoscopic case, however, angular symmetry is usu- 
ally broken by the confining potential. Assuming chaotic 
classical dynamics, the magnitude of the unperturbed wave- 
functions, |0 P | 2 , and the corresponding derivatives, |</jJJ 2 , 
at a given position are statistically independent and both 
obey the Porter- Thomas distribution. 20,24 The perturbed wave- 
functions and their derivatives can then be found from the 
transformation Eq. (|5). Thus, in a mesoscopic situation, the 
MND response does not vanish even at the K-edge. We shall 
see below that this indeed leads to qualitative differences in 
comparison with the bulk behavior: We predict a slightly 
peaked, rather than rounded, K-edge in generic nanosystems. 
At the L-edge, there is a strong MND response, similar to that 
in the bulk metallic case, because the dipole matrix element is 
directly proportional to the amplitude of the perturbed wave- 
function at the position of the perturbation. 
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FIG. 2: (Color online) Chaotic single-particle wave-functions sub- 
ject to a rank one perturbation, (a) Unperturbed and (b), (c) per- 
turbed wave-function probabilities for i = 0, 1, 2 and M, i.e., for the 
lowest three eigenstates and the state at the Fermi energy Ef- We 
assume a two-dimensional chaotic system (TV = 100, M = 50) with 
unperturbed energies at their mean (bulk-like) values and model the 
unperturbed wave-functions as random superpositions of 100 plane 
waves. The intensities along a line that contains the perturbation, 
located at f c = 0, are shown. The normalization volume of the wave- 
function is Q — N. In (b), a weak perturbation causes only slight 
changes in the intensities. In contrast, a strong perturbation in (c) 
causes the wave-function intensity jV'ol 2 corresponding to the bound 
state to pile up at the position of the perturbation. Screening of the 
core hole is done by the bound state on a length scale of order the 
Fermi wavelength, \f, indicated by the black bar. 




III. BOUND STATE 

When the perturbation strength exceeds a certain value, ap- 
proximately v c /d < —1, the lowest perturbed energy level 
Ao will have an energy significantly below all the other lev- 
els (Fig. [T). As discussed in Ref. 5 the average position Ao is 
given in this regime by 

Nd 

A °- e ° = - e d/M_l ■ (10) 
and its fluctuations are negligible. 



As discussed by Friedel 25 and others ^ 26 ' 27 the existence of 
this low energy level is associated with the formation of a 
bound state, which completely screens the core hole pertur- 
bation potential. To illustrate this we employ the Berry- Voros 
conjecture 28 ^ to model the unperturbed wave-functions of 
the chaotic mesoscopic system as a random superposition of 
plane waves. Figure[2fa) shows wave-function intensities for 
both the low-lying eigenstates and those at the Fermi level. 
The final wave function intensities are then shown upon ap- 
plication of a weak [Fig. |2|b)] or strong [Fig. |2|c)] rank one 
perturbation. Whereas the small perturbation has little effect 
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on the wave-function probabilities, the situation changes strik- 
ingly when a strong perturbation is applied: the lowest per- 
turbed eigenstate ipo concentrates at the position of the pertur- 
bation. The extension of the peak in IV'ol 2 is of order the Fermi 
wavelength \p, which for systems with electronic densities 
corresponding to r s > 1, to which our rank one perturbation 
model applies, is of order the screening length A scrccn . 

The existence of a bound state has several important con- 
sequences. The first is the existence of a secondary band of 
absorption corresponding to final states for which the bound 
state is unoccupied (formation of a shake-up pair involving the 
bound state). The secondary band is well separated in energy 
from the main band (where the bound state is filled in the final 
state), and we shall not consider it here in much detail. 

Another consequence of direct relevance to our study is the 
effect of the bound state on the magnitude of the dipole matrix 
elements because of the large amplitude ijjo{r c )- For strong 
perturbation, the probability density |-0o(r c )| 2 at r c is signif- 
icantly larger than the sum for all the other states (we return 
to this point below). In contrast, this is not the case for the 
derivative of ipo at r c , since the largeness of the amplitude 
is compensated by the fact that ipo possesses a maximum in 
the vicinity of f c . Therefore, even for strong interaction, the 
bound state is not going to play a dominant role for the K-edge 
absorption spectra. On the other hand, for the L-edge, where 
the dipole matrix element is proportional to ip Kt replacement 
processes through the bound state will, for strong perturba- 
tion, dominate the absorption. 

In light of this discussion of the bound state, it is instructive 
to see how the Friedel sum rule is treated in our model Hamil- 
tonian approach, Eqs. (f2]l-((5]), in the semiconductor case in 
which there is a single type of conduction electron. Indeed, 
in this model, the Coulomb interaction between electrons is 
neglected, but the fact that the charge of the core hole is +e 
is taken into account by requiring that the phase shift at the 
Fermi energy caused by V c is So = —ir/2. In this way each 
spin channel provides half a charge to screen the core hole. 

Looking at the energy dependence of the phase shift for our 
model, one realizes, however, that the manner of this screen- 
ing is somewhat unintuitive. The most natural initial supposi- 
tion is that the phase shift decreases smoothly from at the 
bottom of the band (ie. little effect of the perturbation on 
states at the band edge) to — tt/2 at the Fermi energy. How- 
ever, what happens in practice is that the phase shift starts 
at — 7r at the bottom of the band and increases gradually to 
—tt/2 at the Fermi energy. In other words, the bound state 
provides a charge — 2e because both spin species have to be 
taken into account, which means that the core hole is actu- 
ally over-screened. All the other (perturbed) wave-functions, 
instead of participating in the screening, are actually pushed 
away from the location of the perturbation, providing an ef- 
fective charge +e near the core hole, leaving, of course, the 
required net screening charge of — e. 

A natural question is whether the scenario described above 
actually happens in physical (experimental) systems, where 
Coulomb interaction between the electrons and chemistry at 
the lattice level both occur. It is clear, for instance, that the 
two electrons occupying the bound state interact strongly with 



each other. If this interaction energy is larger than the dif- 
ference between the bound state and Fermi energies, it will 
prevent double occupation of the bound state and give rise to 
a local moment and so Kondo physics. For the experimen- 
tal realizations we have in mind, however, the bound state 
wave-function will be spread out on the scale of the Fermi 
wavelength (it is not a deep level of an impurity), and simple 
estimates show that a local moment is not expected. 

Assuming Kondo physics is not involved, we now have to 
consider how "real" the bound state actually is in practice and 
if its properties are the same in experimental systems as in 
our model. This question is, of course, not specific to the 
mesoscopic problem. As early as 1952, Friedel discussed 
in detail the physical reality of the bound state in rank one 
models 25 (see also the discussion by Combescot and Nozieres 
in Ref.26|). For example, the screening of the core hole at 
the lithium K-edge is done by the 2s conduction electrons, 
while the sodium L^^-edge is screened by its 3s conduction 
electrons^ It is precisely those s-orbitals which, according to 
the picture developed here following Ref. 25, take the role of 
a bound state. 

Thus the boundstate is a physical reality. But one should 
bear in mind that its extension in space, which controls the 
size of the dipole matrix element uj c o, can be heavily influ- 
enced by the local chemistry or other factors not described by 
our model. In the two limiting cases - (1) if the absorption 
is entirely dominated by the bound state (L-shell with strong 
perturbation), or (2) if the bound state is playing a negligible 
role (weak perturbation or K-shell) - the fact that oj c q may 
not have the proper physical value is of little relevance: In the 
former case, only an overall prefactor is involved in which we 
are in any case not interested, and in the latter case, an incor- 
rect magnitude will obviously not affect at all the description. 
Note furthermore that in the strong perturbation limit, since 
the phase shift at the Fermi energy is essentially independent 
of v c for v c /d <C —l,v c can be chosen so that the energy of 
the bound state is the physical one. 

The situation is more complicated for intermediate situa- 
tions, where the processes involving the bound state have a 
similar contribution to the absorption spectra as those not in- 
volving the bound state (L-shell with v c /d ~ —1). Since our 
model has only one parameter, it cannot reproduce arbitrary 
values of both the phase shift at the Fermi energy and the ra- 
tio between cj c o and the mean value of the other dipole matrix 
elements. In such a situation, it is then necessary to look into 
the details of the bound state for the physical system before 
employing our model. We shall come back to this point in 
Section IVIIII when discussing particular mesoscopic realiza- 
tions. However, in the following, we mainly discuss the limit- 
ing cases for which the physical relevance is unambiguous. 



IV. METHOD 

We will see in this section that the absorption amplitude 
A(uS) for the L-edge is determined entirely by the unper- 
turbed and perturbed eigenvalues {e} and {A}. For the K- 
edge, knowledge of the derivative of the unperturbed wave- 
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functions at r c is also required. For a chaotic system, these 
latter are statistically independent of the spectra and obey a 
Porter-Thomas distribution, 24 and thus do not pose any partic- 
ular difficulty. To study numerically the (statistical) proper- 
ties of A(u)), we therefore generate one particle spectra with 
the known joint distribution 23 by using a Metropolis algorithm 
and deduce from them all the statistical quantities needed to 
characterize the absorption spectra. (Some details about the 
algorithm we have used are given in footnote 30,) This Sec- 
tion derives the basic expressions needed for this purpose. 

Photo-absorption cross section. — Following the Fermi 
golden rule based treatment by Tanabe and Ohtaka, 3 we write 
the photo-absorption cross section A{u>) in terms of the ma- 
trix element of the dipole operator D between the unperturbed 
many -body ground state with the core level c filled, and 
energy E§, and the perturbed (final) many-body state with an 
additional conduction electron \ f) at energy Ef = Eq+Hu, 



2tt 



/ 



2 6{E f - E^-huj) . (11) 



For clarity, we consider spinless electrons for now, and re- 
turn at the end of this Section to the modification introduced 
by spins. The unperturbed ground-state, therefore, comprises 
M electrons on levels to M — 1 with the core level filled, 



1*8) = 



M-l 

n« 



10} 



(12) 



the contribution of the direct process can be expressed as 
A d {uj th ) = A°(w th )|A| 2 . This makes evident the role of AOC 
and the fact that A d vanishes in the thermodynamic limit. 

Even at threshold the direct term is not the only contribution 
to the absorption amplitude: since c K ^ c K , the terms k < M 
in Eq. (fT~4T > are non-zero. These are known as replacement 
processes, Fig. [Tib). In terms of the generalized overlap, 



A, 



(*o|ctc 7 |$ ) 



(16) 



of the unperturbed ground-state with the perturbed state in 
which the particle in the orbital p, < M — 1 has been pro- 
moted to the orbital 7 > M, the total photo-absorption cross 
section at threshold reads 



A d - r (u; th ) 



2tt 

T 
Y 



M M 

°| n ^ e w ^i c c n c i 4 

^=0 |U=0 i=M-l 



W cM L 



M-l 



a i-E 



M=0 



w cM A 



(17) 



Direct and replacement contribution away from 
threshold. — For higher photon energies, such that the 
final state is obtained by adding a particle in the orbital 
7 > M to the perturbed ground state *S?q, the l ast equation is 
readily generalized by the substitution M — > 7, 



The core electron is created (annihilated) by the operator 
cj(c c ). In the perturbed final states \^ f) there are M+l con- 
duction electrons and the core level is empty, 



i*/> = n 4io) 



k filled 



The dipole operator in second quantized form is 



iV-l 



D = ^ w CK (4 Cc + h.c. 



{=0 



(14) 



A d > r (uj) 



2tt 



w C7 A| 1 



M-l 



(18) 



(23) The replacement overlap enters the photo-absorption cross 
section via the ratio w CAl A^ 7 /(w; C7 A) in the second term. 
This ratio can be expressed as a product of eigenvalue differ- 
ences using Eqs. (15)-(19) of Ref.|5j. To this end, we need the 
dipole matrix elements w cfl discussed in Sec. Ill B I Recalling 
that w CM = ip'p (f c ) = u'^ at the K-edge, and = ip^ (f c ) = 
at the L-edge, we express the ratio in Eq. (TT~8b as 



the dipole matrix elements w CK were discussed in Sec. Ill B I 

Photo-absorption processes at threshold. - Let us begin 
our discussion with the absorption at the threshold energy ui t h- 
Then, the only possible final state is the perturbed ground state 
with the core electron excited to level M just above the Fermi 
energy; no shake-up processes are possible. The direct pro- 
cess, Fig. [Tf a), is defined by keeping only k = M in the dipole 
operator dT4b - the term which acts between the core electron 
and the lowest unfilled level. The contribution of the direct 
process is, then, 



(15) 



|C C1 . . . C M |w cM C M C c| c M-l 



444io) 



A d (io th ) 
~h~ 

= ^ Km I 2 1 (Wo) I 2 oc Km| 2 |A| 2 . 

Introducing for reference the amplitude for the bare process 
in which many-body effects are ignored, A°(u} t h) oc |w c m| 2 , 



^c/iA^ 7 zi^Ap 7 

— j u. 



M-l 



w CJ A 



where 



11 w 7 A 



K=0 



A K — A 7 
A K — A,, 



M-l 

n 

fc=0 



A, 



£fc — A 



1 



K-edge 
L-edge 



7 
(19) 

(20) 



carries the symmetry-dependence of the dipole matrix ele- 
ments. The structure of Eq. ( TT8l suggests furthermore the 
introduction of a function p(j), 



(21) 



which can be computed using (\9l . 
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Shake-up processes. — When the energy of the incident 
photon is at least two level spacings above the threshold 
energy, particle-hole pairs can be created in the final state, 
Fig. \V[c). in the one-pair shake-up contribution, two elec- 
trons are excited above Ep in levels 71 and 72 with the level 
fi2 empty in the final state. This situation can be handled in 
close analogy to the replacement process above by introduc- 
ing a renumbered level sequence: let {A*} be a renumbering 
of { A K } in which the level A M2 is skipped whereas A 7l and A 72 
are appended as elements X\ J _ 1 and A M . Let {u 1 } and {u 11 } 
be similarly renumbered sequences. Then, the the one-pair 
shake-up photo-absorption cross section is^ 



Fig- Ed). Note that the energy of the incident photon can be 
shared between the two spin channels. 

In practice, to account for the spin of the electrons the 
photo-absorption spectrum for the optically active channel has 
to be convoluted with curves for the spectator channel. This 
will lead to a slight smearing-out of features obtained for the 
photo-absorption cross section for the optically active chan- 
nel, as we will see below. 



THE SIGNIFICANT REPLACEMENT AND SHAKE-UP 
PROCESSES 



A sh (iu) = 

where 

p 1 (7 2 ) = w cj2 



2tt 



E 

{^2,71,72}!,, 



^(72) 



(22) 



M-l 

E 

fj.=0 



M-l 



n n 



\\, T-r £fc - K 



k=0(/m) 



W TT 

A* 11 e k -X\ { 



f k=Q 



(23) 

Here, the factor T 1 generalizes Eq. ( f2Qb and takes the values 



T 1 



Umu'^/uIu'Ij 



K-edge 
L-edge . 



(24) 



Note that A sh (uj) does not change when 71 and 72 are inter- 
changed in the above equations. 

Generalization to cases with two and more shake-up pairs 
is straightforward: In Eq. ( 1221 ). the overlap of the initial state 
with a state with two or more shake-up pairs is needed. The 
corresponding functions p 2 , p 3 ,... are obtained based on 
renumbering the energy levels such that an index shift occurs 
for each empty level below Ep, and the filled levels above are 
appended. 

As the energy above the threshold increases, the number of 
energetically allowed final states with an arbitrary number of 
particle-hole excitations increases exponentially, and their ex- 
haustive enumeration quickly becomes a hopeless task. How- 
ever, as we shall demonstrate in the next section, the number 
of final states that actually contribute to the absorption pro- 
cess remains finite, and actually not very large. This is what 
makes our approach practical in the end. 

Spin, the spectator channel. — We end this Section by dis- 
cussing how the above picture is modified when the spin of 
the electrons is taken into account. Let us choose the axis of 
spin quantization such that the core electron excited into the 
conduction band has spin up. Since the dipole operator D 
is spin independent, all the discussion above concerning di- 
rect, replacement and shake-up terms applies to the excited 
spin channel. The electrons with spin opposite to the excited 
spin - referred to as the spectator channel - are not connected 
by D to the core electrons; however, they are affected by the 
core hole potential, and their energies and wave-functions are 
modified. As a consequence, the ground state is subject to the 
Anderson orthogonality catastrophe, and some excited states 
may have non-zero overlap with the unperturbed ground state, 



We have seen that it is straight forward to express the ma- 
trix elements appearing in the Fermi golden rule approach, 
Eq. ( ITTb . in terms of the one particle energies {e} and {A}. 
However, the total number of final states increases exponen- 
tially with the energy above threshold. It is therefore neces- 
sary to identify more precisely which final states do actually 
contribute to absorption, and show that the number of such 
states is not prohibitive. 

Toward this end, Fig. [3] shows the spectral weight of the 
unperturbed (many-body) ground state in the perturbed basis, 
as a function of the perturbed state's energy. In addition, the 
cumulative spectral weight (summation up to the perturbed 
state energy) is shown. Including all possible terms and all 
energies, the cumulative spectral weight will, of course, be 
|$o| 2 = 1- Fig. [3] tells us that in practice not all but rather 
only a few terms are needed to reach a spectral weight of 1. 
In particular, it is not necessary to include shake-up processes 
of all orders and energies: Fig.[3]suggests that including terms 
with up to three shake-up pairs and energies up to one and 
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FIG. 3: (Color online) Spectral weight of the unperturbed many-body 
ground state j$o) in terms of perturbed many-body states 
classified by their (average) energy from threshold measured in units 
of mean level spacings d. (v c /d = -10, M — 50,N = 100, COE 
statistics.) The threshold for the secondary band (right panel), when 
the bound state is empty in the final state, is Md greater than the 
threshold with bound state filled (left). The lower (thick blue) curve 
is the energy-resolved width of |$o) in the perturbed basis, equiv- 
alent to the contribution of the spectator channel to the absorption 
cross section. The upper (thinner red) curves show the cumulative 
spectral weight taking into account terms with one, up to two, and up 
to three shake-up pairs (dashed, dotted with symbols, and full lines, 
respectively). Remarkably, less than 0.1% of the weight is missed 
when including only up to two shake-up pairs. The slow saturation 
of the total weight (taking place on the scale of the band width) is a 
characteristic of AOC. 
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a half the band width is 100% sufficient and including terms 
with up to two shake-up pairs captures more than 99.9% of 
the weight. Indeed, one-pair shake-up processes provide the 
dominant contribution to the spectral weight 31 (about 92% for 
v c = —10). The calculation of the photo-absorption spectra 
below includes processes involving up to two shake-up pairs. 

The thicker blue curves in Fig. [3] show the energy distribu- 
tion of the spectral weight (processes with up to three shake- 
up pairs are included). Evidently, a significant percentage of 
the total weight is borne by states of low energy (near the exci- 
tation threshold), or by states with energy equal to or slightly 
above the energy Et> s — Aj\/+i — Ao necessary to promote the 
bound sate electron into an empty orbital. 

Because the dipole operator D is a one-particle operator, a 
non-negligible matrix element (>I'/|Z)c ( t| < I > o) requires \^ f) = 
c^l^y) where n can be arbitrary but \^) is restricted to be 
one of the states which overlaps significantly with <£>o. There- 
fore, the total number of final states that need to be considered 
grows only quadratically with the energy above threshold, and 
there is no problem in enumerating them all. 
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FIG. 4: (Color online) One-pair shake-up and replacement overlaps 
for the bulk-like case, (a) Intermediate perturbation strength, v c /d = 
—0.5. (b) Strong perturbation, v c /d=— 10. For almost all (/i, 7) the 
replacement overlap |A^ 7 | 2 is zero. Non-zero values arise for (1) 
replacement through the bound state with the excited electron close 
to Ef (p, = and 7 > Af, black arrow), and (2) shake-up pairs 
formed in the vicinity of the Fermi edge (/1, 7 both close to M, red 
arrow). For a strong perturbation, replacement through the bound 
state becomes the dominant process. 



For processes involving only one shake-up pair, a more de- 
tailed representation of the decomposition of the unperturbed 
ground state is shown in Fig.|U where |Ap 7 | 2 is shown as a 
function of p and 7 for the bulk-like case. These bulk-like 
overlaps, which also determine the replacement contribution, 
provide a good estimate for the mean value in the mesoscopic 
case and are useful to estimate the relative importance of the 
different processes. We see immediately that most of these 
overlaps are very small. There are two notable areas of excep- 
tion: one where the particle-hole pair lies close to the Fermi 
energy (/1, 7«M), and a second for terms involving the bound 
state. For a small perturbation, the Fermi energy peak dom- 
inates. However, as soon as a bound state develops upon in- 
creasing the perturbation, the overlaps |A^| 2 start to grow 
and eventually overwhelm those near the Fermi energy. 



VI. PHOTO-ABSORPTION SPECTRA: AVERAGED 
CROSS SECTION 

In the mesoscopic case, the photo-absorption threshold en- 
ergy fluctuates from sample to sample. We assume here that 
Wth can be determined experimentally, and that energies and 
spectra are then measured with respect to this energy. There- 
fore, we will often give the photon excess energy relative to 
the threshold energy in mean level spacings, fi (uj — uo t h)/d. 

Furthermore, in a mesoscopic system the final state en- 
ergy Ef can only take discrete values. A mesoscopic photo- 
absorption spectrum is therefore comprised of a series of 5- 
peaks (broadened in experiment). We are interested in par- 
ticular in the prefactor given by the Fermi golden rule matrix 
elements. First, we discuss their average in this Section, turn- 
ing to fluctuations in the next. In all cases, we present only the 
photo-absorption from the primary band, for which the bound 
state is occupied in the final state. The absorption is normal- 
ized to the total absorption of the primary band. 



A. K-Edge 

Fig.[5]shows the average photo-absorption cross section for 
a K-edge. The time-reversal non-invariant case (CUE) is con- 
sidered [see Ref. 4jfor a similar illustration of the time-reversal 
invariant case (COE)]. Note that the replacement contribution 
decreases rapidly away from threshold, at which point shake- 
up processes become important. The large replacement ampli- 
tude at threshold causes the mesoscopic K-edge to be peaked. 

Fig-Hb) shows the modification made by spin. The con- 
tribution of electrons with spin opposite to the excited core 
electron - the spectator channel - is shown. The full spin 
cross section is obtained after convolution of the spectra of 
the two spin species. The slight peak at the K-edge threshold 
is maintained in the full spin spectrum. 

Comparison of the mesoscopic and bulk-like situations is 
shown in Fig. [6] At a K-edge (lower curves), the mesoscopic 
and bulk-like photo-absorption spectra are qualitatively dif- 
ferent. The bulk-like K-edge shows the rounded behavior 
expected from AOC, though the threshold value is non-zero 



9 




<co - co th ) / d (go - co th ) / d 



FIG. 5: (Color online) Average mesoscopic absorption spectra at the 
K-edge (N = 100, M = 50, v c /d = -10, CUE), (a) The total ab- 
sorption cross section in the active channel (full line) is the sum of 
the direct/replacement (triangles) and shake-up (diamonds) contribu- 
tions. For comparison, the bare contribution (dashed-dotted line) and 
the direct process alone (squares) are shown, (b) Active (full circles) 
and spectator (open circles) channel spectra separately, as well as the 
full spin photo-absorption cross section obtained after convolution 
(down-triangles). The edge is slightly peaked. 



FIG. 7: (Color online) Average mesoscopic absorption spectra at the 
L-edge (N= 100, M = 50, v c /d = - 10, CUE), (a) The spectrum of 
the optically active channel is the sum of the direct/replacement (tri- 
angles) and shake-up (diamonds) contributions, (b) Active (full cir- 
cles) and spectator (open circles) channel spectra separately, as well 
as the full spin photo-absorption cross section obtained after convo- 
lution (down- triangles). The contribution of the spectator channel is 
identical to that in Fig. Ob). The peak at the L-edge is much more 
pronounced than that at the K-edge (Fig.O and extends over several 
mean level spacings in photon energy. 
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FIG. 6: (Color online) Average mesoscopic (triangles) and bulk-like 
(quadrangles) spectra as a function of energy from threshold at both 
a K- and L-edge (TV = 100, M = 50, v c /d = -10, CUE). Whereas 
the bulk-like and mesoscopic-chaotic results coincide for the L-edge, 
there is a clear difference in the K-edge spectra: the bulk-like edge 
is rounded whereas a generic mesoscopic system yields a slightly 
peaked edge on average. The dashed curves are the mesoscopic spec- 
tra in a COE situation; they are nearly indistinguishable from the 
CUE case. 



due to incompleteness of AOC in a finite system. The aver- 
age mesoscopic K-edge spectra, on the other hand, is slightly 
peaked at threshold. However, as the photon energy becomes 
only a few mean level spacings above threshold, the average 



mesoscopic K-edge photo-absorption quickly approaches the 
bulk-like result from above. 

The reason for the different K-edge behavior in bulk-like 
and mesoscopic samples lies in the dipole matrix elements. In 
the bulk-like situation, dipole selection rules cause the matrix 
elements in the I = channel to vanish, and therefore there is 
zero MND response. In contrast, the dipole matrix elements 
in the mesoscopic situation are generally non-zero random 
numbers with a Porter-Thomas distribution, originating from 
the distribution of wave-function derivatives as discussed in 
Sec. Ill B I Therefore, there is a MND response near thresh- 
old in the mesoscopic situation that counteracts the AOC edge 
rounding, thus causing the peaked K-edge. 



B. L-Edge 



The dipole matrix elements w cfl for the L-edge are propor- 
tional to the wave-function ip fJ ,(F c )=U tl at the position of the 
perturbation (see Sec. Ill Bl ). We therefore find that whenever 
the perturbation is strong enough to form the bound state this 
latter will play a very significant role. Indeed, in the situation 
where the / = channel provides all the screening of the im- 
purity (i.e. <5o = 7r/2), the photo-absorption process is entirely 
dominate dby the term uj c qcIc c in the dipole operator Eq. ( TBl i. 



10 




band filling (co - co, h +0.5d) / N d band filling (co - co th +0.5d) / N d 

FIG. 8: (Color online) Mesoscopic averaged photo-absorption cross section as a function of the number 2M of particles in a half-filled band 
for (a)-(b) K-edge and (c)-(d) L-edge, and both weak coupling [(a),(c) v c /d= — 0.3] and strong coupling [(b), (d) v c /d= — 10]. Results are for 
the COE with full spin up to excitation energies a quarter of the band width above threshold, and are normalized by the bare photo-absorption 
value. The K-edge appears, apart from the behavior directly at threshold, rounded, and the rounding increases for more particles in the system. 
In this sense AOC wins the competition with MND as the thermodynamic limit is approached. However, the slight peak at the edge persists 
as the signature characteristic of a mesoscopic-coherent system. The L-edge clearly is peaked, and this peak sharpens with increasing particle 
number. For strong coupling, the L-edge is completely dominated by replacement through the bound state, making the magnitude at the L-edge 
much larger. 



The shape of the average photo-absorption spectra can in this 
case be essentially understood from the energy dependence of 
the overlap Aq 7 shown in Fig. [4] a sharp peak for 7 = M is 
followed by a relatively long tail. As shown in Figs. [7] this 
is the behavior of the L-edge photo-absorption spectrum in 
the strong perturbation regime. Note that the L-edge peak is 
considerably stronger than for the K-edge (compare to Fig.[7j>. 
Convolution of the sharp peak in the active channel with the 
spectator spin result still yields a prominent peak. 

In comparing with the bulk-like case, we start by empha- 
sizing again that the physics for an L-edge in the strong per- 
turbation regime is dominated by the bound state. As there 
are no strong differences between the bound state in the bulk- 
like and mesoscopic cases, we expect the results to be similar. 
Fig. [6] shows, in fact, a stronger result: for a strong perturba- 
tion, the mesoscopic and bulk-like spectra at an L-edge agree 
quantitatively^ 



C. Dependence on Perturbation Strength and Number of 
Particles 

So far our major focus has been the strong perturbation 
regime and a model nanosystem with 100 electrons (50 elec- 
trons per spin species in a half-filled band). Since the pro- 
cesses that determine the shape of the edge, namely the AOC 
and MND responses, depend on the number of particles, we 
will now address how the (average) photo-absorption spectra 



change as the number of particles in the system is varied. This 
is of particular interest with respect to experiments since the 
number of electrons in the system can be adjusted by means 
of the geometry (size), gate voltages, or the density of states 
(doping). It is convenient at the same time to vary the strength 
of the perturbation produced by the core hole. The strong per- 
turbation regime describes semiconductor heterostructures. In 
these systems, only s-conduction electrons are present and 
available for screening; the Friedel sum rule then implies the 
strong perturbation regime. In contrast, Fermi edge singu- 
larities in metals are described by weaker perturbationsr^ 
corresponding to a smaller phase shift for the s-electrons and 
in agreement with the fact that other channels of conduction 
electrons (p, d, ...) are involved in the screening. Metallic 
nanoparticles are one mesoscopic system where similar val- 
ues for the phase shifts occur. Other situations where the small 
perturbation regime is relevant could arise. 

Fig. [8] shows photo-absorption spectra at both the K- and 
L-edge for two different perturbation strengths with the num- 
ber of electrons ranging from 24 to 200. (Note that COE 
statistics are used here rather than the CUE statistics used in 
Figs. 017]) The weaker perturbation (v c = —0.3) produces a 
phase shift at Ep which is typical of a metallic environment, 
while the larger strength produces complete s-wave screening 
(<5 f ~ — 7r/2) suitable for semiconductors. All curves show 
the main FES signatures: for a K-edge, a peak at threshold 
superposed on a rounded edge, while at an L-edge, a strong 
peak at threshold. Clearly, the FES signatures are enhanced 
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when the perturbation is strong. 

For the weaker perturbation, the bound state is not formed. 
Absorption near the edge is nevertheless enhanced because of 
correlation between the spectrum and the values of the wave- 
functions at r c . This implies in particular that all the terms 
in the replacement sum Eq. ( fTTl ) have the same sign, as in the 
bulk and despite the random character of the wave-functions. 

Since the mean level spacing d = 2n/N is TV-dependent, 
the energy from threshold is given here in units of band filling 
above threshold, (uj — ci>th) / N d (such that an excitation into 
the highest level 7 = N— 1 corresponds to a value 1 /2) in order 
to allow for a direct comparison of the curves. Note that we do 
not confine our attention to the immediate threshold vicinity 
but rather consider excitation energies that allow electrons to 
fill states up to 3/4 of the band width. 

We first discuss the K-edge spectra in Fig. [8] First, note that 
all four curves converge at large energy to the same behavior. 
In fact, they approach a line corresponding to the bare process 
(~ |w C7 | 2 ). This is easily understood: Once the photon en- 
ergy is well above threshold, so that u> — ui t h is significantly 
larger than the width of the unperturbed ground state in the 
perturbed basis (see Fig. [3), the final states are necessarily of 
the form \^ f) = cj t |\f r °) where \^) is one of the states with 
a significant overlap with $o- All the "extra energy" is borne 
by the highly excited one particle state k. In this case, only 
the term w CK c£.c c of the dipole operator D contributes to the 
absorption - there are no replacement processes. Summing 
over final states amounts to averaging lu^ k on an energy win- 
dow equal to the width of $ m the perturbed basis, and so the 
result is the same as for the bare process. 

Moving toward the K-edge threshold, one sees that the 
photo-absorption is first suppressed and then jumps right at 
threshold. The suppression is simply the manifestation of the 
AOC familiar from the bulk K-edge. In fact, as N increases, 
the rounding becomes more pronounced. When the perturba- 
tion is weak, the resulting points lie right on the bulk power- 
law curve. Thus, for weak perturbation, the average meso- 
scopic photo-absorption in the N — > 00 limit yields the bulk 
singularity. For strong perturbation (the semiconductor case), 
there are substantial deviations from the bulk power-law: at 
large energies the average mesoscopic photo-absorption is 
suppressed while for energies just above threshold, as for the 
threshold point itself, photo-absorption is enhanced. 

A striking difference between the bulk and mesoscopic K- 
edge spectra occurs right at threshold: rather than following a 
rounded edge, the average mesoscopic response shows a peak 
(see also Fig. |6). The magnitude of this peak does depend 
on the strength of the perturbation but is approximately N- 
independent. For the strong perturbation case (v c = — 10), the 
peak is approximately 50% larger than the second point and 
about 3 times larger than the bulk result. This enhancement 
near the Fermi edge comes from the dipole matrix elements 
which, for the generic nanosystems considered here, are non- 
zero even at the K-edge: breaking of rotational symmetry al- 
lows MND processes in the mesoscopic case while only AOC 
is present in bulk. 

The L-edge spectra in Fig. [8] show the strongly peaked edge 
characteristic of the MND singularity. In the case of a weak 



perturbation, panel (c), there is extremely good agreement be- 
tween the average mesoscopic photo-absorption and the bulk 
power-law response for all energies and all N. Right at thresh- 
old, it appears that the average mesoscopic result is slightly 
larger than the bulk-like result. The peak becomes more pro- 
nounced as N increases because, as for the K-edge, one is 
able to access smaller energies with respect to the band width. 

For strong perturbation at the L-edge, note the very large 
magnitude of the photo-absorption. This stems from the 
fact that replacement processes through the bound state com- 
pletely dominate. In this case, the bulk power law provides a 
qualitative guide to the peak, but, as for the K-edge, does not 
agree quantitatively. 

VII. FLUCTUATIONS OF THE PHOTO-ABSORPTION 
CROSS SECTION 

Turning from the average photo-absorption, we now inves- 
tigate a quantity inherent to mesoscopic systems and a key 
characteristic of the photo-absorption cross section: its fluctu- 
ations. Fig.|9]shows several examples of the photo-absorption 
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FIG. 9: (Color online) Individual K-edge absorption spectra of 
six mesoscopic samples (marked by different symbols), illustrating 
the outcome expected in real single-sample measurements (iV = 
100, M = 50, v c /d=— 10, COE). Fluctuations occur in both the en- 
ergy and the cross section. The spectra consist, strictly speaking, of a 
series of 5-functions (experimentally broadened) whose weight and 
positions are marked by the points; the lines are guides for the eye. 
(a) Photo-absorption cross-section as a function of photon energy to. 
The cross section is highly fluctuating, not necessarily peaked, and 
can even zigzag. Both the threshold energy and the distance between 
the points vary on the scale of the mean level spacing, (b) Spectra 
shifted such that their threshold energies coincide. A peaked edge 
visibly emerges. 
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FIG. 10: (Color online) Fluctuations of the photo-absorption cross 
section at the K-edge for (a) strong and (b) weak perturbation (N = 
100, M — 50, CUE, optically active spin). The average photo- 
absorption cross section as a function of energy, (A(ui)), is shown 
for processes with different average excitation excitation (marked 
ojth + rf, <^th + 2d, . . .). Note that energies are measured with respect 
to the threshold energy wth (vertical line) as would be the case in 
experiments. The behavior of the threshold energy and cross section 
are shown in the inset: both are approximately Gaussian distributed. 
The peak next to threshold is clearly asymmetric with the maximum 
photo-absorption cross section found at energies distinctly below 
the average value. The curves broaden and symmetrize away from 
threshold. The area under the curves is the total photo-absorption; 
the (slight) peak at threshold is seen by noticing that the area in the 
inset is largest, as well as being greater for larger perturbations. 

cross-section of individual systems. As for the average, all 
data presented are for absorption in the primary band and 
are normalized to the total absorption of the primary band. 
There are two sources of fluctuations: wave-function (matrix- 
element) fluctuations and energy-level fluctuations. 

First, the wave-function amplitudes at the location of the 
core hole vary, causing the dipole matrix elements to fluctu- 
ate. This is the most dramatic mesoscopic effect, causing large 
fluctuations in the photo-absorption cross section (Pig. [9j». In 
particular, the shape of the spectrum of an individual system is 
not necessarily peaked, but can be "rounded", almost uniform, 
or even zigzag. 

Second, the energy level effect creates fluctuations in the 
photo-absorption as well as in the photon energies that can be 
absorbed. The former is small compared to the effect of the 
wave-functions, but the latter can be quite significant. The 
absorbed photon energies fluctuate by an amount of order the 
a mean level spacing d, with increasing width further away 
from threshold. Because core electrons at different locations 
can be excited by the different photons, subsequent measure- 
ments of one and the same system will result in different final 
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FIG. 11: (Color online) Fluctuations of the photo-absorption cross 
section at the L-edge for (a) strong and (b) weak perturbation (N = 
100, M = 50, CUE, optically active spin). Explanations are the same 
as for Fig. [10] The strongly peaked FES is evident. 

(perturbed) energy levels {A}, and therefore give rise to dif- 
ferent spectra even though the unperturbed energy levels {e} 
remain the same. With single measurements of different sys- 
tems, as we assume here, both the {e} and {A} would vary. 

Note that the threshold energy varies from dot to dot, re- 
sulting in a relative shift of the spectra of a few mean level 
spacings [Fig. [9ja)]. In Fig.[9jb) the spectra are shown with 
respect to their threshold energy, clearly enhancing the organi- 
zation of the "measurement" points. Still, the fluctuations add 
up away from threshold such that the spectra can span rather 
different energy ranges. 

In order to analyze the fluctuations more systematically, 
Figs. [10] and [TTJ show the photo-absorption cross section as 
a weighted histogram, highlighting the fluctuations in en- 
ergy weighted by absorption intensity. The total cross sec- 
tion is broken down by grouping processes which have the 
same mean photon energy, and then for each of these groups 
the average (A(ui)) is plotted. Two perturbation strengths 
are contrasted, and energies are measured with respect to the 
threshold energy (as would be the case in experiments). Note 
that the curves for the different classes have substantial over- 
lap. For the peak next to threshold, higher cross sections are 
achieved at smaller energies for all perturbation strengths and 
at both the K- and L-edges. Further away from the Fermi 
edge, the histograms become symmetric, broader, and more 
Gaussian, as one expects from the Central Limit Theorem. 

In the rest of this Section, we focus on the statistical prop- 
erties of the absorption amplitude. We aggregate all the pro- 
cesses (final states) taking place at the same mean energy; the 
resulting probability distributions are shown in Fig. [12] 

Most of the curves closely resemble the Porter- Thomas dis- 
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FIG. 12: (Color online) Distribution of the mesoscopic photo- 
absorption fluctuations (iV = 100, M = 50, CUE): for a strong per- 
turbation v c /d — —10, (a) the K-edge and (b) the L-edge, and for a 
weak perturbation v c /d= — 0.3 at (c) the L-edge with the K-edge in 
the inset. The probability distribution of the photo-absorption A(ui) 
normalized by its mean value (A(ui)) is shown for different mean 
excitation energies (tj), both near threshold (curves with symbols) 
and further away from threshold (full line) where a Gaussian shape 
emerges. For comparison, a Porter-Thomas distribution is indicated 
by the dashed line. Near threshold, the distributions are Porter- 
Thomas in almost all cases because the fluctuations in the dipole ma- 
trix elements dominate and cancel all correlations from the overlap. 
An exception occurs at threshold for the L-edge where the excited 
electron sits in the first level above the Fermi energy. In this case, the 
distribution resembles that of a ground state overlap. 

tribution, in particular for the weaker perturbation (exceptions 
are discussed below). This originates from the dipole matrix 
elements ui CK being proportional to Porter-Thomas distributed 
wave-function derivatives (K-edge) or wave-functions (L- 
edge) at the position of the perturbation (Sec. Ill BY The cor- 
responding randomness overwhelms any correlations in the 
ground state, replacement, or shake-up overlaps. Further away 
from threshold, as more and more processes are included, the 
distribution becomes more Gaussian. 

At the L-edge in the limit of strong perturbations, we find 
however clear deviations from the Porter- Thomas distribution 
[Fig. [T2l" b)l. At threshold, (u) — (u)th)> the corresponding 



curve (circles) is very similar to the probability distribution 
P(| A| 2 ) of the ground state overlap (between the unperturbed 
and perturbed many-body ground state) presented in the first 
paper of this series (see Fig. 3 of Ref.|5|). The origin of this 
behavior can be traced back to the fact that, as discussed 
in Sec. [Hi] the absorption in the L-edge strong perturbation 
regime is entirely dominated by processes involving the bound 
state. Therefore at threshold 

A(iu th )~ Ko| 2 |A 0M | 2 . (25) 

As |if c o| 2 does not fluctuate, the distribution of A(uj t \i) is the 
same as that of |A uM | 2 , which is equivalent 5 to an Ander- 
son ground state overlap for positive perturbation and phase 
shift 7T— \Sf\ [Eq. (21) in Ref. 5]. In the strong perturbation 
limit, \6f\ =7r/2 = n— \Sp\, and therefore the distribution 
of |Aqm| 2 is the same as the ground state overlap distribu- 
tion. For weaker perturbations, such that the bound state still 
dominates the absorption at threshold but with a phase shift 
\Sf\ < 7r/2, the symmetry relation between the Anderson 
overlap for negative perturbation and the replacement over- 
lap involving the bound state for a positive perturbation still 
holds, but the negative and positive phase are no longer equal. 

The mechanism underlying the fluctuations in the L-edge 
strong perturbation case when the photon energy is above 
threshold deserve further discussion. On the one hand, as il- 
lustrated in Fig. [3] the overlap between final states f =c\^® 
and the unperturbed ground state rapidly becomes extremely 
small. On the other hand, the bound state dipole matrix ele- 
ment w c o is much larger than the others (basically |u> c o| 2 is of 
the same order of magnitude as l^cil 2 )- It is therefore 

not clear a priori which effect, the smallness of the overlaps 
or the largeness of the matrix element, will dominate. 

The answer to this question is given in Fig. Qj] which shows 
the relative importance of bound state dominated processes in 
various cases. We see, first, that for weak perturbations, the 
importance of the lowest perturbed one-particle state (which 
is not here properly speaking a bound state) is, as expected, 
marginal. We further see that for the K-edge with strong per- 
turbation, this contribution is not negligible (about 25% near 
the edge for the case considered) but neither is it dominant. 
Furthermore, it decays rapidly to zero as the photon energy 
becomes larger than the width of the unperturbed ground state 
in the perturbed basis. 

For the strong perturbation L-edge case, panel (d), which is 
our main interest here, we see that the largeness of the dipole 
matrix element is the dominant effect, even at large photon 
energy. The overlap does not decrease fast enough to com- 
pensate for this, and only a moderate decrease of the relative 
contribution of the boundstate is observed. As a consequence, 
there are no fluctuations associated with the dipole matrix ele- 
ment (since |iu c rj | does not fluctuate) — the fluctuations derive 
entirely from those of the overlap between the unperturbed 
ground state |$o) an d the final states \^ /) from which the 
bound state has been removed. As the energy of the photon in- 
creases, the number of shake up processes (i.e. of final states) 
increases, and the fluctuations become more Gaussian. 

So far we have discussed fluctuations in the CUE situa- 
tion where time-reversal symmetry is broken, corresponding 
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FIG. 13: (Color online) Contribution of bound state processes to the 
average photo-absorption cross section as a function of photon en- 
ergy above threshold (N = 100, M = 50, CUE, optically active 
spin). All replacement and shake-up processes involving the bound 
state are included. Upper (lower) panels: weak (strong) perturba- 
tion. Left panels: K-edge. Right panels: L-edge. For weak perturba- 
tion, the contribution is, as expected, small. For strong perturbation, 
processes through the bound state make a significant but decaying 
contribution at the K-edge; at the L-edge, they dominate. 



to the average photo-absorption spectra shown in Sec. [VI] 
Very similar results apply to the COE case (presence of time- 
reversal symmetry) as we found above for the average photo- 
absorption. The major difference comes from the change in 
the Porter-Thomas distribution (it now diverges at zero); once 
this is taken into account, the fluctuation characteristics (to the 
extent described here) are the same for COE and CUE. 



VIII. EXPERIMENTAL REALIZATIONS 

We end this paper by discussing some possible experimen- 
tal realizations of the physics we have described. Beyond a 
simple illustrative purpose, this Section provides an opportu- 
nity to discuss in more detail the applicability of the model we 
have used to real physical systems. 

The model in Eqs. ©-(O has been built to contain the main 
ingredient of the physics under consideration, namely the sud- 
den appearance of a local perturbation whose phase shift at the 
Fermi energy is determined by the Friedel sum rule. It turns 
out, however, that the output of the model does not depend 
only on its properties at the Fermi energy, even qualitatively. 
In particular, the bound state created in the strong perturbation 
regime plays a major role in influencing both the mean and 
fluctuations of the photo-absorption spectrum. It is therefore 
useful to discuss for a few examples how much this mecha- 
nism is expected to be taken into account correctly. 

Since the basic physics we describe - a localized (rank one) 
perturbation acting on a Fermi sea of chaotic electrons - is 




FIG. 14: (Color online) Experimental arrangement to test our predic- 
tion of a transition from a rounded K-edge to a slightly peaked edge 
as the system size is diminished to the mesoscopic coherent regime. 
We propose a quantum dot array in a semiconductor heterostructure 
where the plane of the 2DEG also contains impurities that provide a 
localized state in the GaAs band gap. 



very general, it can be realized in various systems: 

(1) An obvious experiment would be the direct realization 
of x-ray absorption measurements in small metallic nanopar- 
ticles. Either colloidal nanoparticles with diameter 1-2 nm or 
"metallic molecules" such as AU55 might be used. The techni- 
cal requirements for such an experiment are beyond standard 
capability at present: Although the spatial resolution of x-ray 
micro-beams has been demonstrated at low temperature, their 
energy resolution seems to be insufficient. 33 Such an experi- 
ment may, of course, be possible in the future. 

With regard to the importance of the bound state (discus- 
sion at the end of SectionlHlli. the bound state should play the 
same general role in a nanoparticle as in the bulk material. 
In the bulk metals that have been studied most carefully rill 
several angular momentum channels are involved in screen- 
ing, yielding phase shifts less than 1 . These, then, are in the 
weak perturbation regime: the bound state plays a minor role 
in both the bulk and mesoscopic situations. 

(2) Another, intrinsically different, example occurs in a 
double quantum dot where the tunneling of an electron into 
(or out of) the system causes a rank one perturbation due to 
the constriction which mediates a sudden change of the wave- 
function in one dot to the other^i^ The AOC analysis in 
Ref. |34| shows that the phase shifts involved are ±7r/4. Such 
a phase shift is produced by v c /d in the range —0.3 to —0.5. 
Thus for this situation as well, the perturbation is weak, and 
the effect of the bound state is minor. 

(3) Finally, photo-absorption in quantum dots is a third ex- 
ample, one which is feasible with existing standard semicon- 
ductor heterostructure technology. 36 Since for semiconductor 
conduction electrons all of the screening comes from a single 
channel, the phase shift is ir/2, and the role of the bound state 
is critical. We now describe in detail this possibility. 

First, a quantum dot array can be formed in a semiconduc- 
tor heterostructure like GaAs/AlGaAs using standard growth 
and etching techniques (see Fig. [14}. The dot diameter can 
be made smaller than 100 nanometers, with a variation of 
±5 nanometers, using state-of-the-art techniques. This gives 
a mean level spacing of about 30 ± 3 /ieV. Slight fabrication 
imperfections guarantee that the dots are chaotic, all slightly 
different but with a well defined average. The layer of the two- 
dimensional electron gas (2DEG) cuts through all the dots 
which "stand" as pillars on a GaAs substrate. Each dot con- 
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tains several tens of conduction electrons; this number can be 
varied in the fabrication process by adjusting the doping. 

The core electron of the original x-ray edge problem is now 
replaced by an electron in an impurity level within the GaAs 
bandgap. In order to realize the assumptions made in the pa- 
per, the electron has to be in a localized state, and in order to 
probe a K-edge situation, it has to have s-like symmetry. This 
can be realized using a suitable deep level, or by doping with 
N (isoelectronic to As). The impurities host a s-type bound 
state occupied by an electron that is subsequently excited into 
the conduction band. The variety of available dopants should 
even allow one to change the symmetry of the impurity level, 
and thus to study the L-edge behavior. The impurities have to 
be placed into the 2DEG in order to ensure a sufficient overlap 
and have a suitable density (few impurities per dot). 

The excitation of the impurity-pinned electrons is achieved 
using a micrometer laser (rather than an x-ray). The resolution 
of these lasers is around 10 /ieV, thus below the mean level 
spacing of the electrons in the quantum dots and sufficient to 
resolve the expected effects within a few mean level spacings 
from threshold. The photo-absorption cross section has to be 
measured at low temperature. 

The first quantity investigated in this paper is the aver- 
age photo-absorption cross section. Although it is tempting 
to measure it directly as the photo-absorption of the array, 
this may not be possible (yet) due to fabrication uncertainty. 
The Fermi energy in each dot (measured from the impurity 
level) will vary by about a mean level spacing from dot to 
dot due to single particle effects; in addition, the charging en- 
ergy may lead to a slightly larger variation, perhaps two mean 
level spacings. Finally, the mean level spacing itself will vary 
by about 10% from dot to dot (assuming 5% uncertainty in 
the linear dimension). Taken together, the dot-to-dot varia- 
tion may wash-out the expected threshold peak in the photo- 
absorption, which extends, after all, over only about two mean 
level spacings. Therefore, each dot may have to be measured 
individually. Afterwards, the average over an ensemble (the 
array) has to be performed. In this way, information about the 
fluctuations can, of course, be extracted as well. 

IX. CONCLUSIONS 

In this paper we have studied Fermi edge singularities in 
the mesoscopic regime, in particular photo-absorption spectra 
of generic, chaotic-coherent mesoscopic systems with a finite 
number of electrons. The basic underlying physics that we 
study is, however, much more general, and the model can be 
applied to any situation where a sudden, localized (rank one) 
perturbation acts on a finite number of chaotic electrons via a 
(dipole) matrix element. 

The photo-absorption cross section differs from the naive 
expectation due to two counteracting many-body responses 
resulting from the sudden perturbation of the system when 
a core or impurity electron is removed: the Anderson or- 
thogonality catastrophe (AOC) and the Mahan-Nozieres- 
DeDominicis (MND) contribution. We have studied AOC in 
detail in the first paper of this series. 5 We found, for instance, 



that typical mesoscopic systems are rather far from the ther- 
modynamic (bulk) limit and AOC is incomplete. This is also 
reflected in the photo-absorption spectra: the K-edge photo- 
absorption cross section near threshold remains finite. 

Whereas the threshold behavior in x-ray photoemission (ex- 
cited electron leaves the system) can be deduced from the 
power law exponent of the AOC overlap alone, this does 
not apply to photo-absorption or luminescence experiments 
(where the electron is excited into, or relaxes from, the con- 
duction band) where, besides AOC, the MND response is cru- 
cial. It counteracts the AOC effect if the dipole selection rules 
are fulfilled. In the present paper we included this response 
and computed photo-absorption spectra of generic meso- 
scopic systems in various situations using a Fermi golden rule 
approach. In metals the result is typically a rounded K- and a 
peaked L-edge in the photo-absorption cross section. 

We confirm the peaked L-edge behavior in the average 
photo-absorption cross section for mesoscopic systems. At 
the K-edge, however, we find characteristic changes that al- 
low one to infer the presence of chaotic-coherent dynamics 
of electrons in nanosystems from the near-threshold behavior 
of the average photo-absorption. As the system size is dimin- 
ished to reach the mesoscopic-coherent scale, the most impor- 
tant change occurs in the dipole matrix elements. The chaotic 
dynamics of the (conduction) electrons makes them distinc- 
tively different from the bulk-like situation: The electronic 
wave-function and its derivative at a certain position are inde- 
pendent, and the dipole matrix elements at the K-edge (oc ip' K ) 
are non-zero. As a consequence, mesoscopic K-edge spectra 
are peaked and not rounded as in the metallic (bulk-like) case. 
This peak is visible right at threshold and is a direct signature 
of chaotic-coherent dynamics of electrons in mesoscopic sys- 
tems. We propose in detail experiments where this prediction 
can be tested using current semiconductor technology. 

The fluctuations of the photo-absorption cross section are 
largely dominated by the Porter- Thomas distribution of the 
dipole matrix elements. However, deviations occur at the 
L-edge right at threshold when the cross section probability 
distribution resembles that of the AOC ground state overlap. 
The reason for this behavior is the dominance of replacement 
through the bound state, a recurring theme in this subject. 
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